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Abstract. We provide an efficient reduction from the problem of querying approximate multiplicatively 
weighted farthest neighbors in a metric space to the unweighted problem. Combining our techniques 
with core-sets for approximate unweighted farthest neighbors, we show how to find (1 + e) -approximate 
farthest neighbors in time O(logrc) per query in £>-dimensional Euclidean space for any constants D and 8. 
As an application, we find an 0(n\ogn) expected time algorithm for choosing the center of a star topology 
network connecting a given set of points, so as to approximately minimize the maximum dilation between 
any pair of points. 

1 Introduction 

Data structures for proximity problems such as finding closest or farthest neighbors or maintaining 
closest or farthest pairs in sets of points have been a central topic in computational geometry for a 
long time [2]. Due to the difficulty of solving these problems exactly in high dimensions, there has 
been much work on approximate versions of these problems, in which we seek neighbors whose 
distance is within a (1 +e) factor of optimal [4]. In this paper, we consider the version of this 
problem in which we seek to approximately answer farthest neighbor queries for point sets with 
multiplicative weights. That is, we have a set of points pu each with a weight w(pi), and for any 
query point q we seek to approximate max/ w {pi)d(pi, q). 
We provide the following new results: 

- We describe in Theorem 3 a general reduction from the approximate weighted farthest neighbor 
query problem to the approximate unweighted farthest neighbor query problem. In any metric 
space, suppose that there exists a data structure that can answer unweighted (1 +e) -approximate 
farthest neighbor queries for a given «-item point set in query time Q(n,e), space S(n,e), and 
preprocessing time P(n,e). Then our reduction provides a data structure for answering (1 + 
e) -approximate weighted farthest neighbor queries in time 0(logn + (Mog \)Q(n,z/2)) per 
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query, that uses space 0(n + (-log-)S(n,E/2)), and can be constructed in preprocessing time 
0(nlogn+(Ilog±)P(n,e/2)). 

- We apply core-sets [1] to find a data structure for the approximate unweighted farthest neighbor 
query problem in R D , for any constant D, with query time 0{z~), space 0(e~), and pre- 
processing time 0{n + zi~ D ) (Theorem 5). Applying our reduction results in a data structure 
for the approximate weighted farthest neighbor query problem in MP, with preprocessing time 
0(rclog7i + £2~ D log^), query time 0(logn + e^~ log^), and space 0(n) (Corollary 6). 

- As a motivating example for our data structures, we consider the problem of finding a star- 
topology network for a set of points in R D , having one of the input points as the hub of the net- 
work, and minimizing the maximum dilation of any network path between any pair of points. By 



results of Eppstein and Wortman [5], this problem can be solved exactly in time 0(n2 a ^ log 2 «) 
in the plane, and 0(n 2 ) in any higher dimension. By using our data structure for approximate 
weighted farthest neighbor queries, we find in Corollary 16 a solution to this problem having di- 
lation within a (1 + e) factor of optimal, for any constant £ and constant dimension, in expected 
time 0(nlogn). More generally, as shown in Theorem 15, we can approximately evaluate the 
dilations that would be achieved by using every input point as the hub of a star topology, in 
expected 0(n\ogn) total time. 

2 Problem Definition 

We first define the weighted farthest neighbor query problem in a metric space, and then extend 
its definition to relevant variants and special cases. Let M = (S,d) be a metric space. For any two 
points p and q in S, d(p,q) denotes the distance between them. We assume we are given as input a 
finite set P C S; we denote by the points in P. In addition, we are given positive weights on each 
point of P, which we denote by a function w : P i— ► M + . We wish to preprocess the point set such 
that for any query point q £ 5, we can quickly produce a point / £ P that maximizes the weighted 
distance from q to r. More precisely, / should be the point in P that maximizes (or approximately 
maximizes) the weighted distance d w (q,pt) = w{pi) -d(q,pj) from the query point q. If we restrict 
the weights of all points to be equal, then this problem simplifies to the unweighted farthest neighbor 
query problem. 

Finding the exact farthest neighbor to a query point can be computationally expensive. Hence, 
we are interested in approximate solutions to this problem, with sufficient guarantees that our al- 
gorithm produces a reasonably distant point in the weighted sense. To achieve this, we define the 
approximate weighted farthest neighbor query problem, wherein we are given an input of the same 
type as the weighted farthest neighbor query problem, and wish to preprocess the point set so that 
we can efficiently produce an £-approximate weighted farthest point. That is, we wish to find a point 
r <EP such that d w (q,r) > (1 —e)d w (q,pi) for some predefined £ and all p\ £ P. A similar definition 
holds for the approximate unweighted farthest neighbor query problem as well. When the metric 
space under consideration is Euclidean, we call the problem the Euclidean approximate weighted 
farthest neighbor query problem. 

Without loss of generality, we can assume that < w{p) < 1 for all p £ P, and that exactly 
one point p\ has weight 1. The first assumption can be made to hold by dividing all weights by the 
largest weight; this leaves approximations to the weighted farthest neighbor unchanged. The second 
assumption can be made to hold by arbitrarily breaking ties among contenders for the maximum- 
weight point. 

3 Reduction from Weighted to Unweighted 

First, let us suppose that we have a family of algorithms JZ e for the approximate unweighted farthest 
neighbor query problem defined over any metric space M = (S,d). A e consists of two components, 
one for preprocessing and the other for querying. For an input instance of the approximate un- 
weighted farthest neighbor query problem, let / £ P maximize the distance between our query point 
q and /. Using A e , we can find an £-approximate farthest point r such that (1 — e)d(q,f) < d(q, r). 
The running time of A e has two components, namely preprocessing time and query time denoted 
P(n,e), and Q(n,e) respectively. It has a space requirement of S(n,e). Our real concern is the 
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weighted version of the problem and in this section, we provide a data structure to reduce an in- 
stance of the approximate weighted farthest neighbor query problem such that it can be solved by 
invoking A E . We call this family of reduction algorithms 2^. 

The preprocessing and querying components of %£ are provided in Procedures 1 and 2 respec- 
tively. For convenience, we separate out the points in S into two subsets Si and S2. The first subset Si 
contains all points in P with weights in [e/2, 1], and S2 contains P\S\. We use different approaches 
for the two subsets producing three candidate approximate farthest neighbors. Our final result is the 
farther of the three candidates in the weighted sense. We now show that the algorithm is correct both 
when f £ Si and when / £ S2. 



Case / € Si : 

In this case, we only consider points in Si . Therefore, we can assume that the weights are within 
a factor of 2/e of each other. We partition the points set Si into buckets such that each bucket 
Bj contains all points in Si with weight in [(e/2)(l + e/2)(''~ 1 >,(e/2)(l +e/2)'). The number of 
buckets needed to cover Si is ®((l/e)ln(l/e)). 

Our data structure consists of separate instances (consisting of both the preprocessing and 
querying components) of A e with the error parameter e/2 for each bucket Bj. In addition, each 
Bj is preprocessed by the preprocessing instance of A E . To query our data structure for the (e/2)- 
approximate farthest neighbor kj € Bj from q, we run the querying instance of A e on each prepro- 
cessed bucket Bj. Our candidate approximate weighted farthest neighbor for this case will be the 
n = argmaxfc; <%,&,)• 

Lemma 1 If f G Si, then 

d w (q,n)>(l-e)d w (q,f), (1) 
where f is the exact weighted farthest point and r\ is as defined above. 

Proof. Let us just consider bucket Bj whose weights are between some (e/2)(l +e/2)(''- 1 ) and 
(e/2)(l +e/2)'. Let f be the exact weighted farthest point from q in Bj. Within the bucket Bj, kj is 
an (e/2)-approximate farthest point if we don't consider the weights. I.e., 

(l-e/2)d(q,fj)<d(q,kj)- 

However, when we consider the weighted distance, w(kj) and w{f) differ by a factor of at most 
(1 +e/2) and hence our inequality for the weighted case becomes 

1 - e/2 

dw{qji) < d w (q,kj). 



l+e/2 

_ \ (■{ cA fnr O \ A ,t,c 

l+e/2 



Since j— 4j > (1 — e) for e > 0, we get 



(l-e)d w (cj,f)<d w (q,ki). (2) 

Let the bucket Bj* contain the weighted farthest point /. By definition, d w (q,r\) > d w (q,kj*), 
which proves the lemma when combined with Equation 2. □ 
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Procedure 1 Preprocessing for the approximate unweighted farthest neighbor problem 



1: INPUT: Set of points P with n — \P\ and a distance metric d. A fixed error parameter e. 

{Farthest point is in 5i .} 
2: Let Si CPbe points with weights in [2/e, 1]. 

3: Partition Si into buckets Bq,B\,... such that each B, contains all points in S\ whose weights are in ((e/2)(l + 

£/2r 1 ,(e/2)(l + (£/2))'']. 
4: Preprocess each bucket by calling the preprocessing component of A e instance with error parameter e/2. 

{Farthest point is in 52-} 
5: Assign p\ to be point in P with highest weight (assumed to be 1). 
6: Sort S2 = .P \ Si in non-decreasing order of d(p\,s) for all s € S2. 
7: for each s taken in reverse sorted order from S2 do 
8: if 3x y s in S2 such that d w (p\ ,s) < d w {p\ ,x) then 
9: S 2 <-S 2 \{s}. 

{Note that this can be done in 0(n) time by updating the x that maximizes d w (p\,x) at each iteration.} 
10: end if 
11: end for 

12: Ensure that S2 is stored in a manner suitable for binary searching. 



Procedure 2 Query step of the approximate unweighted farthest neighbor problem 

1: INPUT: Query point q, B, for all 1 < i < l/eln(l/e), and S 2 . 

{Farthest point is in S\ .} 
2: Assign r\ <— q 
3: for each B, do 
4: fci<-Jl( e/ 2)(fli) 
5: if dyv(q,ri) <d w {q : k i ) then 
6: r\ <- kj 
7: end if 
8: end for 

{Farthest point is in S2.} 

9: Binary Search through S2 for point ro such that d(p\ , r^) >d[jp\,q)jz and d w {p\ ,r2) is maximized. 
10: Return r *- arg max {n ripi y (d w (q,n), d w (q,r 2 ), d„ (q,pi)) 
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Case f £ S 2 : 

We sort the points in 52 in non-decreasing order of their distance from p\\ recall that p\ is the only 
point in P with w(pi) = 1. We then pare down 52 by preserving points x such that d w {p\,x) > 
d w (p\,y) for all y >- x in the sorted sequence of points, and discarding the rest. Our data structure 
for this case consists of the remaining reduced sorted list, stored in a manner suitable for performing 
binary searches (e.g. an array or binary search tree). To query the data structure, we binary search 
for the point r2 that is the first point in the reduced sorted S2 farther than d(pi,q)/e from p\. Our 
two candidates for this case are r2 and p\. 

Lemma 2 If f € 52, and p\ is not an e-approximate weighted farthest point from q, then d w (q, r 2 ) > 



Proof. Consider the ball whose center is p\ and radius is 2d(q,p\)/e. Since / € 52, we can say that 
d(qj) >2d{q,pi)/e. Otherwise, d w (q,f) < (e/2)d{q,f) <d(q,pi) =d w (q,p 1 ) (since w(pi) = 1) 
and this contradicts our assumption that / £ 52- In other words, both / and r 2 fall outside this ball, 
the former by argument and the latter by construction. Hence, 



In addition, r 2 is chosen over / as being the farthest weighted neighbor of p\ because (by construc- 
tion) it is the weighted farthest point from p\ outside the ball in consideration. Hence, 



(l-e)d w (q,f). 



d(pi,r 2 ) > 2d(pi,q)/e. 



(3) 



d w (pi,r 2 ) > d w {p u f). 



(4) 



Consider Apiqr 2 . 



d{p\,r 2 ) < d(p u q)+d(q,r 2 ) 

< (e/2)d(p u r 2 )+d{q,r 2 ) 
=> d(q,r 2 )>(l-e/2)d( Pl ,r 2 ). 



(due to 3) 



Again, 



d(p u r 2 )+d(pi,q) > d(q,r 2 ) 
d(p u r 2 ) + (e/2)d(pi,r 2 ) >d(q,r 2 ) 
(\+E/2)d( Pl ,r 2 )>d(q,r 2 ). 



(due to 3) 



Therefore, 



(l-e/2)d( Pl ,r 2 ) <d(q,r 2 ) <(l+£/2)d(p u r 2 ). 



(5) 



With similar treatment of Apiqf, we get 



{\-z/2)d{p u f)<d(qJ) < (l+£/2)d( Pl J). 



(6) 



From Equation 6, we have 



d w {qj) < {\+z/2)d w { Pl J) 
< (l+e/2)d w (p u r 2 ) 



(From 4) 




(From 5) 



=> d w (q,r 2 ) > (l-e)d w (q,f). 



□ 
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Theorem 3 The family of reduction algorithms %£ answer (1 + €)-approximate weighted farthest 
neighbor queries in time f?(log?i + (^log -)Q(n,e/2)) per query, use space 0(n+(-log-)S(n,e/2)), 
and the data structure can be constructed in preprocessing time 0(nlogn + (^log ^)P(n,e/2)). 

Proof. The proof of approximation guarantee follows from Lemmas 1 and 2 in a straightforward 
manner. Preprocessing takes time 0((^log ^)P{n,z/2)) for partitioning into buckets and calling the 
preprocessing component of & e , and 0(n\ogn) for sorting the points, hence accounting for the total 
preprocessing time. While querying, we call the querying component of A E for each bucket costing 
us 0((ilog ^)Q(n,e/2)), and the binary search takes time 0(logn) adding to the total querying 
time. Similarly, the space taken to store the preprocessed buckets and the sorted sequence of points 
is 0((^-log ^)S(n,e/2)) and 0{n) respectively, adding to the total space required. 

□ 

The preprocessing and querying components (Procedures 1 and 2 respectively) are presented in 
pseudocode format. Each procedure has two parts to it corresponding to the two cases when / G S\ 
and f £ S 2 - In Procedure 1, we first partition Si into buckets and preprocess them. Secondly, we 
sort S 2 and pare it down according to the requirements of Lemma 2. In Procedure 2, we choose 
the farthest neighbor r\ of q in S\ by querying each bucket and choosing the farthest from the pool 
of results. Secondly, we binary search in the reduced S 2 to obtain the point r 2 that is at or just 
farther than a distance of 2d(q,pi )/z from q. Our final £-approximate farthest neighbor r of q is the 
weighted farthest point in {r\ , r 2 , p\ } from q. 

4 Euclidean Approximate Unweighted Farthest Neighbor Queries 

In this section, our problem definition remains intact except for the restriction that S is M°, where 
D is a constant, with the Euclidean distance metric. Notice that our points are unweighted and we 
are seeking a family of approximation algorithms j? £ that we assumed to exist as a black box in the 
previous section. Let / € P maximize the distance between our query point q and points in S. For 
a given £, we are asked to find an £-approximate farthest point r such that d(q,r) > (1 —e)d(q,f). 
We can solve this using the £-kernel technique surveyed in [1]. Let u be the unit vector in some 
direction. We denote the directional width of P in direction u by w(u,P) and is given by 

w(u,P) = max(M,s) — mm(u,s), 

SEP SEP 

where (u,s) is the dot product of u and s. An £-kernel is a subset K of P such that for all unit 
directions u, 

(l-e)w(n,P) <w(u,K). 

It is now useful to state a theorem from [1,3] that provides us an £-kernel in time linear in n and 
polynomial in 1 /e with the dimension D appearing in the exponent. 

Theorem 4 Given a set Pofn points in MP and a parameter £ > 0, one can compute an e-kernel of 
P of size 0(\/e^ D - 1 ^ 2 ) in time 0(n + 1 /e^ 3 / 2 )). 

Consider points si and S2 in P and &i and Ic2 in K that maximize the directional widths in the 
direction u, which is the unit vector in the direction of qf. 

(1 -e)w(u,P) = (1 — e)(u,s\) - (1 -e)(u,s 2 ) < (u,ki) - {u,k 2 ). 
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Now, (u,k 2 ) > (u,S2), because otherwise, w(u,P) can be maximized further. Hence with some sub- 
stitution, we get 

(1 -e)(w,Ji) < {u,h). 

The point / maximizes the left hand side in the direction of qf because it is the farthest point 
from q. Therefore, the above inequality suggests that k\ € K is an £-approximate farthest neighbor 
of q. This implies that the point farthest from q obtained by sequentially searching the £-kernel K 
will be an £-approximate farthest neighbor of q. Hence, we can state the following theorem. 

Theorem 5 Given a set P of n points in M° and a parameter £ > 0, there exists a family of algo- 
rithms A z that answers the e-approximate unweighted farthest neighbor query problem in prepro- 
cessing time 0{n + e^ ~ D ), and the query time and space requirements are both 0{z~ ). 

Combining this theorem with our general reduction from the weighted to the unweighted prob- 
lem gives us the following corollary. 

Corollary 6 Given a set P of n points in M. D , for any constant D, and given a weight function w : 
P i — > M + and a parameter £ > 0, there exists a family of algorithms that answers the e-approximate 
weighted farthest neighbor query problem in preprocessing time 0(nlogn + £2~ D log -) and query 
time 0(log n + £T~ log i ) with a space requirement of 0(n). 

Proof. We apply Theorem 3 to reduce the weighted problem to a collection of instances of the 
unweighted problem, which we solve using Theorem 5. The query time comes from combining 
the results of these two theorems, and the space bound comes from the fact that our overall data 
structure consists only of a sequence of points together with a family of disjoint subsets of points 
forming one core-set for each of the instances of the unweighted problem used by Theorem 3. To 
calculate the time bound, note that Theorem 3 takes time O(nlogn) to sort the low-weight points by 
distance from p\ , and in the same time we may also perform the partition of high- weight points into 
buckets used by that theorem. Once we have partitioned the points into buckets, we construct for 
each bucket an instance of the data structure of Theorem 5; this takes time 0{n\ + e^~ D ) per bucket 
and adding these bounds over all buckets gives the stated preprocessing time bound. □ 

5 Constrained Minimum Dilation Stars 

The dilation between two vertices v and w of a weighted graph is defined as the ratio of the weight 
of the shortest path from v to w, divided by the direct distance between v and w. The dilation of the 
entire graph is defined as the greatest dilation between any pair of vertices in the graph. A star is a 
connected graph with exactly one internal vertex, called its center. Any collection of n points admits 
n possible stars, since any individual point may serve as a star center. In this section we consider 
the problem of computing the dilation of all of these n stars. A solution to this problem provides the 
foundation for a solution to the problem of choosing an optimal center: simply search for the point 
whose corresponding dilation value is smallest. 

Eppstein and Wortman considered [5] the problem of selecting star centers that are optimal 
with respect to dilation. They showed that for any set of n points P C MP, there exists a set C of 
0(n) pairs of points such that the worst pair for any center c G U° is contained in C. They give an 
0[n log n) -time algorithm for constructing C, and go on to consider the problem of computing the 
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center that admits a star with minimal dilation. They present an 0(nlogn) expected-time algorithm 
for the case when c may be any point in MP and D is an arbitrary constant, and an 0(n2 a ^ log 2 72) 
expected-time algorithm for the case when c is constrained to be one of the input points and D = 2. 

These results imply that the dilation of all n stars with centers from P may be computed in 0(n 2 ) 
time: construct C, then for each c G P and (v,w) G C, evaluate the dilation of the path (v,c,w). In 
this section we improve on this time bound through approximation. We will show that a (1 — e)- 
approximation of the dilation of all n stars may be computed in 0{n\ogn) expected time, and hence 
an approximately optimal center c G P may be identified in 0(n\ogn) expected time. 

Our approximation algorithm first uses the results of [5] to compute the optimal Copt G ^P- 
Note that Copt may or may not be in P, and in general will not be. Our algorithm then partitions 
P into two subsets: the k G 0(1) points nearest to Copt, and the other n — k points. We say that 
a star center c G P is k-low if it is one of the k points nearest to Copt, or k-high otherwise. The 
dilation values for all the k-low centers are computed exactly in 0(nlogn) time using a combination 
of known techniques. 

The dilation between any two points Pi,pj G P through any &-high center c may be approximated 
by a weighted distance from c to the centroid of p\ and pj. Hence the the dilation of the star centered 
on any /c-high c may be approximated by the distance from c to the centroid farthest from c. We use 
the data structure described in Section 4 to answer each of these n — k weighted farthest-neighbor 
queries in O(logrc) time, which makes the overall running time of the approximation algorithm 
0(n\ogn). 

Definition 7 Let G be a Euclidean star with vertices V and center c. The dilation between any 
v,w G V \ {c} is 



8(c,v,w) 

and the dilation of the star is 



\vc\ + \wc\ 



\vw\ 



A(c) = max 8(c,v,w). 

v,weV\{c} 



Definition 8 Let 



— P = (po , . . . ,p n -i ) be a set ofn input points from MP for some constant D, each with the potential 
to be the center of a star with vertices P, 

— c opt G MP be the point minimizing A(copt)> 

— S = (sq,. . . ,s„_i) be the sequence formed by sorting P by distance from Copt, 

— £ > be a constant parameter, 

— r = 2/e-l, 

— k be a constant depending only on T, 

— L = {sj G S I < i < k} be the k-low centers, 

— and H = P\L be the k-high centers. 

We require the following claim, which is proved in [5] : 

Claim 9 Let c be the center of a Euclidean star in MP for D G 0(1) having vertices V. If A(c) < T 
for some constant T, then there exists a constant pr such that for any integer i, the D-dimensional 
annulus centered on c with inner radius p' r and outer radius p'^ 1 contains only 0(1) points from V. 
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Fig. 1. Planar annuli containing the points .v,, v, and w. 



Lemma 10 For any T there exists a constant k depending only on T such that any center Sj G S with 
i > k has dilation A(j,) > T. 

Proof. We first consider the case when A(copt) > T. By definition A(sj) > A(copt) for any Si G S, 
so we have that every A(si) > T regardless of the value of k. 

We now turn to the case when A(copt) < T- Define pr as in Claim 9, and let 

Aj = {x G M. d I p J r < \xc OP t \ < P } r +l } 

be the jth annulus centered on Copt- Define 1 = 1 + [log Pr (r — 1 )] , and suppose Si G Aj. Let v, w be 
two input points that lie in the annulus Aj-i. By the definition of / we have 



2(pf- l - P r / K r 

Since v G Aj-i, we have \sjv\ > (p } r 1 — p J r '), and similarly |j,-w| > (p J r 1 — p J r '). Further, 
\vw\ < 2p/r'. So by substitution 

(M) + (M) >r 



(|vw|) 

A(j f ) > r. 

Thus the lemma holds if the annuli A 7 -,Ay_i, . . . contain at least k points. We can ensure 
this by selecting k such that the points . . . ,Sj necessarily span at least / + 1 annuli. By Claim 9, 
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there exists some m G 0(1) such that any annulus contains no more than m input points. So we set 
k> (1 + l)m. The observation that / depends only on F completes the proof. 

□ 

Corollary 11 For any c Si, A(c) < T; for any c €H, A(c) > T; and \L\ £ 0(1). 

Proof. Each claim follows from Lemma 10 and Definition 8. □ 

Lemma 12 The set L, as well as the quantity A(c)for every c £L, may be computed in 0{n\ogn) 
expected time. 

Proof. The unconstrained center copt may be found in 0(n\ogn) expected time [5]. L may be 
constructed by sorting P by distance from copj, and retaining the first k elements, where k is the 
constant defined in Lemma 10. 

The set C may also be constructed in 0{n\ogn) time [5]. Then any A(c) may be computed by 
evaluating 

A(c) = max 8(c,v,w), 

(y,w)eC 

which takes 0(n) time since |C| € 0(n) [5]. Evaluating all \L\ = k dilation values this way takes 
0(kn) time. So the expected amount of time needed to compute every A(c) is 0(n\ogn+kn), which 
is 0(nlogn) by Corollary 11. □ 

Lemma 13 If c € H, pt,Pj € P be the pair of points such that 8(c,pi,pj) = A(c), and \pic\ > \pjc\, 
then \pjc\ > (1 — e)\pjc\. 

Proof. By the assumption c G H and Corollary 11, 

a, x \p,c\ + \pjc\ 

\PiPj\ 



by the triangle inequality 



\PiPj\ + \Pf\ > \Pi c \ 

\PiPj\ > \pic\-\pjc\ 



\Pic\ + \pjc\ >r 



(\pic\-\pjc\) 

\pic\ + \pjc\ >T(\pic\-\pjc\) 
\pf\+T\pjc\ > T\pic\ - \pic\ 

M(r+i)>| A -c|(r-i) 
, r- 1 , 



By the definition of Y, 
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Fig. 2. Example points p,-, pj, qjj, and i 



. (2/e— 1) — 1 , 
^ * (2/e-l) + l M 



2/e-2 
>(l-e)|pyc|. 



□ 

Lemma 14 Define qtj = ^(pi+pj) andwij = \ ^ p .y then for any c£// and corresponding Pi,Pj € 
P such that 8(c,pi,pj) = A(c), 

Wij-\qijc\ > (l-e)-S(c,pi,pj). 

Proof. By Lemma 13 

\ Pj c\ > (l-e)\ Pi c\ 

and by construction \qijc\ > \pjc\, so 

Wijc\ > (l-e)\pic\ 
2-\qijc\ > (l-e)(|p/c| + |p,-c|). 

By assumption \pic\ > \pjc\, so 

2-\qijc\ > (l-e)(\pic\ + \ Pj c\) 

AW>(1-B)W 
IPiPil PiPj 

and by substitution 

»•/./• 'n.jc > (l-e)-S(c,pi,pj). 

□ 
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Theorem 15 A set ofn values A(p,) subject to 



(l-e)A(p,-)<A(pi)<A(pi) 
may be computed in 0{n\ogn) expected time. 

Proof. By Lemma 12 is is possible to separate the &-low and &-high centers, and compute the exact 
dilation of all the &-low centers, in 0(n\ogn) expected time. Section 4 describes a data structure that 
can answer approximate weighted farthest neighbor queries within a factor of (1 — e) in 0(log?i) 
time after 0(nlogn) preprocessing. As shown in Lemma 14, the dilation function 8 may be approx- 
imated for the fc-high centers up to a factor of (1 — e) using weights wtj = So for each &-high 
Pi, A(pf) may be approximated by the result of a weighted farthest neighbor query. Each of these 
\H\ =n — k queries takes 0(logn) time, for a total expected running time of 0(n\ogn). □ 

Corollary 16 Let Aopt = min p ,ep A(pi). Then a point c G P satisfying 

(1 - e)A pt < A(c) < A 0PT 

may be identified in 0(nlogn) expected time. 

Proof. Theorem 15 shows that the values A(/?;) may be computed in 0{n\ogn) expected time. A 
suitable c may be found by generating these n values, then searching for the smallest A{pi) and 
returning the corresponding pi. □ 
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